New Physics and Neutrino Oscillation 
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Abstract 

Description of neutrino oscillation in the case of Non-Standard neutrino Interaction 
(NSI) is briefly presented. The NSI causes the entanglement between internal degrees 
of freedom of neutrinos (mass, spin, flavour) and other accompanying particles in the 
production and detection processes. In such case neutrinos are mostly in the mixed 
states. Role of the density matrix in description of neutrino oscillation process is shortly 



p i' explained. 



Theory of neutrino oscillation found in late sixties pLj works well and is almost commonly 
accepted. Now days two methods describing particles oscillation are considered, one based 



I on the quantum mechanics (QM) and the second which use the quantum field theory (QFT) 
^ ' (for a review see p]). For the standard relativistic neutrino interaction both methods are 
■ equivalent, and give the same result for final neutrino oscillation rate [3j. In the case of 
^ ■ a NSI there is no, up to now, uniform way of describing the neutrino oscillation process. 

In the standard approach the state of produced neutrinos does not depend on the state of 
accompanying particles, states are disentangled. In such a case neutrinos oscillate in the 
way independent of the spins arrangement of the other particles present in the production 
and detection process. In the QFT approach, and in most studies of the neutrino oscillation 
beyond the Standard Model (SM) [1], oscillation process is assumed to be universal, where 
probability of neutrino transition is possible to be defined and depends only on the initial and 
final states. 

However, beyond the SM the production and detection neutrino states may differ from 
states which are used to construct the Lagrangian of the model before the spontaneous sym- 
metry breaking, the so called flavour states |5]. It can happen also that the right-handed 
neutrino fields appear, although there are a fairly strong limits on couplings involving these 
fields |6j. This causes that neutrino states in the production, detection and propagation pro- 
cesses need not to be given by the simple combination of neutrino mass states with the same 
helicity, 

i i 

which is the basic assumption of QM description of the standard neutrino oscillation. In most 
approaches beyond the SM, instead of the states ([1]) the new one, separate for production (p) 
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and detection (d) are defined 



\<) = j:km)^ wt) = EuiM)^ (2) 
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where the mixing matrices t/^ j and U^^ are connected with the NSI Hamiltonians for pro- 
duction and detection 

m/^\{iyfJp\H''\L;ip)\', \Ulfcx\{h-fD\H^\u,-tD)\'. (3) 

The ip(D) and fp(D) are initial and final states for production (detection) process. Then the 
amplitude of neutrino transition from the initial state \i) to the final state |/) (given by Eqs. 
([I]) or (I2])) is defined in the standard way 

Af, = {f\e-^^'\t), (4) 

where H is Hamiltonian of free neutrinos (in vacuum) or describing neutrino interaction with 
matter particles (in matter). In this description of the oscillation phenomena there is no 
possibility to take into account any kind of intrinsic correlation between neutrinos and other 
accompanying particles, but such internal correlations exist in some models of NSI. 

The natural description of the oscillation phenomena is the QFT approach, where the pro- 
duction and detection neutrino states are not used, all calculations are done in the neutrino 
mass base and it is not necessary to assume that oscillations are described by production 
and detection independent probabilities [7]. In this description there is however one not 
very natural assumption that neutrinos propagate as virtual particles over macroscopic or 
even astronomical distance. Therefore another way of describing the oscillations, without 
the propagation of virtual neutrinos, where the QFT is used to construct the neutrino pro- 
duction and detection states only, has been proposed [8]. Taking into account that the QM 
(internal wave packets) and QFT (external wave packets) descriptions of neutrino oscillations 
are equivalent we proposed the generalization of both approaches where any kind of NSI can 
be taken into account [3] . The proposed description can be summarized in the following points: 

(i) If the M/j^"([p], [A]) are the amplitudes for "a" neutrino produced in the mass "z" and 
helicity " A" state in the process: B + l„ + Ui{X), and the M^f{[q][jl]) describe the "/?" 

neutrino with mass "k" and helicity "/x" in the detection process z/fc(//) + C — )■ + D then 
the amplitude for the a — )■ /3 transition after time T and at a distance L is given by [3] 

iAap{X,fi;[p],[X],[q],[-p])= (5) 

0(T) E M^;-{[%[X])M^f{[q],[jl]) j^^Fp^u{E,{p),p)FnAEk{v).P^^^^^ 



A:=l,2,3 



where the [p] and [A] {[g] and \p]\ describe the average momenta and helicities of all other par- 
ticles in the production (^4, B, /„) {detection (C, D, /^)} process. The functions Fp^k{Ek{p),p) 
and Fp,k{Ek{p),p) are constructed from momenta distribution functions for all particles in 
the production and detection processes respectively without neutrinos (for details see 0) and 



2 



6(T) is the Heaviside step function. 

(ii) The probabihty for the overall process - production, oscillation and detection - is 
calculated as the modulus square of the amplitude ([5]) with proper averaging over initial and 
sum over final particles helicities and averaging over the necessary particles momenta in the 
same Lorentz frame (usually the rest frame of a detector) . Practically, it is more convenient to 
calculate the production and the detection processes separately in their own rest frames, and 
in the final oscillation rates to distinguish the initial neutrino state and the final cross-section 
for the detection process. 

(ii(l)) In the rest frame of A the neutrino initial state is calculated from dynamics of the 
production process A ^ B + la + i^i{X) in the well know way. The neutrino density matrix in 
the mass (z) and helicity (A) base |i, A) is equal 

gl,.^^^,iE,e,ip)= (6) 

where the integral dLips is taken over the part of the phase space, without neutrino energy 
(E) and its momentum direction (6*, the Qx^^^y^ describes the polarization density matrix 
of decaying particle (A) and the factor Na normalizes the density matrix, such that Trg = 1. 
It is convenient to normalize separately the initial state and the final cross section, then for 
the plane wave (not wave packet) approach the final oscillation rate is already normalized 
properly. 

(ii(2)) To get the density matrix in the LAB frame, where decaying particles are moving 
(pions, decaying nuclei, muons), the Lorentz transformation, equivalent to the Wigner rota- 

1/2 l/2t 

tion, Qlab = Dylr QcmD-^ has to be performed. For relativistic neutrinos in the helicity 
base, the matrix — 1, then the plab to a very good approximation is equal to (f^'^ [5] 

qUb{e,, e,, cp) = d'^'qI,,{e, e, ^)d'^'^ ^ qI,, ^ ^"(i?, e, ip), (7) 

where the neutrino energy (E) and its spherical angle {9) are expressed in terms of the 
appropriate quantities E^ and 6b in the LAB frame after the Lorentz transformation. 

(ii(3)) As neutrino detector is usually far away from the production place, to all practical 
purpose it is possible to use average density matrix defined in the following way 

where is the solid angle in which neutrinos "see" the detector. The ^^^^q^ is the neutrino 
energy and angular distribution. For large baseline L, practically 'g°'{Eh) — i?"(6'b = 0, tp). 

(iii) The detection cross section is calculated in the standard way. Such a cross section 
describes the neutrino oscillation and detection process in one formula. Neglecting the co- 
herence effect connected with the particles wave packets, the oscillation cross section is given 
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by 



(9) 
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Y: I dLips Mf/(p;, [71]) g?^,.^,^,{E,) e' 




327rs \pi\ 2sc + 1 



The two formulae (EI) and ([9]) are crucial for our analysis. From (I6j), calculating Tr{{g°')'^), 
we see when the neutrino state is pure, or when it have to be considered as mixed. For 
example we can easy find that in the SM, 



which describes the pure quantum mechanical state. However, there are models of NSI in 
which neutrinos with both helicities can be produced. In such models, in the obvious way, the 
neutrino production states are mixed. But even if only the relativistic neutrinos with negative 
helicities are produced the states can be mixed too, it depends on the specific structure of the 
helicity amplitudes (see e.g. [5]). From the Eq. we can easy find when the final oscillation 
rate factorize and when not. 

In conclusion, we have presented the theory of neutrino oscillation applicable in the case of 
any NSI, where the internal degrees of freedom of neutrinos and other accompanying particles 
in a production and/or in a detection process are entangled. 
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